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r*| ■ Abstract 

a: 

O 1 The contributions of subleading shape functions to inclusive decay distributions of B 

mesons are derived from a systematic two-step matching of QCD current correlators 
onto soft-collinear and heavy-quark effective theory. At tree-level, the results can be 
expressed in terms of forward matrix elements of bi-local light-cone operators. Four- 
^ quark operators, which arise at 0(a s ), are included. Our results are in disagreement 

with some previous studies of subleading shape-function effects. A numerical analysis 
of B — > X u l~v decay distributions suggests that power corrections are small, with the 
possible exception of the endpoint region of the charged-lepton energy spectrum. 



1 Introduction 

Inclusive decays of B mesons into final states containing light particles, such as B — > X u l~v 
and B — > X s 7, play a prominent role in the extraction of the Cabbibo-Kobayashi-Maskawa 
(CKM) matrix element \V U b\ and in searches for physics beyond the Standard Model. Necessary 
experimental cuts in the analysis of these processes restrict the hadronic final state to have 
large energy, Ex ~ mg, but only moderate invariant mass, Mx ~ a/^bAqcd- 111 this region 
of phase space, the inclusive rates can be calculated using a twist expansion, which resums 
infinite sets of power corrections into non-perturbative shape functions PjJ 12 E] . 

The impressive performance of the BaBar, Belle, and CLEO experiments demands a con- 
tinuous effort to reduce the theoretical uncertainties affecting extractions of \V u b\ using inclu- 
sive B decays. Recently, significant progress has been made by systematically incorporating 
higher-order perturbative corrections [31010]. A careful estimate of the theoretical uncertainty 
on \Vub\, extracted using a cut on P + = E x — \Px\, finds an error of 7% due to variations 
of the leading shape function [7j, which can be reduced significantly using information on the 
B — > Xgj photon spectrum. With a 5% relative theoretical error on \V u b\, corrections sup- 
pressed by a power of Aqcd / mb are considered the second largest source of uncertainty. The 
present paper is devoted to a more thorough study of power corrections to inclusive S-meson 
decay distributions in the shape-function region, using the two-step matching procedure de- 
veloped in jUISlin]. Subleading shape functions have been investigated first by Bauer, Luke, 
and Mannel [Sj, and their effects on various inclusive spectra have been analyzed by several 
groups US UH U2] • We disagree with the findings of [HI HO] in some important aspects. 

The hadronic physics in the inclusive semileptonic decay B — > X u l~v is encoded in a 
hadronic tensor W^ u defined via the discontinuity of the forward S-meson matrix element of 
a correlator of two flavor-changing weak currents J M = M7 M (1 — 75) b. We define 



where v is the I?-meson velocity and q the momentum carried by the lepton pair. To be slightly 
more general, we will consider QCD currents of the type jj = bTiq and Jj = qTjb, where 
q is a massless quark, and r$, 1^- are arbitrary Dirac matrices. The corresponding current 
correlator Tij and hadronic tensor Wij can then also be applied to study the contribution of 
the dipole operator Q 77 to B — > X s ^ decay. 

The current correlator receives contributions from different energy scales, which can be 
disentangled using effective field theories. Besides the hard scale m;, and the hadronic scale 
Aqcd, an intermediate "hard-collinear" scale 'nib Aqcd set by the typical invariant mass of 
the hadronic final state is of relevance. Because quantum fluctuations associated with the hard 
and hard-collinear scales can be treated in perturbation theory, the hadronic tensor trivially 
factorizes into products of hard functions (/i ~ m&), jet functions (/i ~ a/^Aqcd) , and 
shape functions defined in terms of 5-meson matrix elements of quark-gluon operators. In 
a first step, hard fluctuations are integrated out by matching the QCD currents onto soft- 
collinear effective theory (SCET) |T51 114j. The current correlator is then expanded in terms 
of light-cone operators in heavy-quark effective theory (HQET) [H)], thereby integrating out 
fluctuations at the hard-collinear scale. The fact that the relevant HQET operators live on 
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the light cone follows from the structure of the multipole expansion of soft fields in SCET. In 
the present work, we carry out the matching procedure at tree level and to order Aqcd/? 7 ^ 
in the heavy-quark expansion. We indicate how our results would change if loop corrections 
were included. 



2 Short-distance expansion of the hadronic tensor 



We begin by recalling some elements of SCET, referring to [M] for a more detailed discussion 
(see also |13j). We work in a reference frame where the B meson is at rest, t> M = (1,0,0,0), 
and where the lepton 3-momentum q points in the negative ^-direction. We introduce two 
light-like vectors n M = (1, 0, 0, 1) and n M = (1, 0, 0, —1), with n ■ n = 2, n ■ v = n ■ v = 1, and 
n = 2v — n. Any 4- vector a M can be decomposed as 

a M = n ■ a + n ' a ~2 a ± = a - + a + + a ± • (2) 

In this basis, v± — and q± = by construction. In the kinematic region of interest, the 
hadronic jet has a hard-collinear momentum, which scales like (P_,P + ,Pj_) ~ m&(l, A, \/A), 
where A ~ Aqcd/^- The momenta of the soft, light constituents of the B meson scale like 
P^ s ~ A QCD ~ m b \. 

SCET is the appropriate effective field theory for the description of the interactions among 
soft and hard-collinear degrees of freedom. Its Lagrangian is organized in an expansion in 
powers of y/\. The leading-order Lagrangian is 

^scet = f 77 f in ' D hc + gn ■ + ip±hc ._ 1 r) ip±.hc J £ 

Z \ in- L>hc J 

+ q ip s q + hiv ■ D s h + £yM > (3) 

where £ is a hard-collinear quark field, q is a soft, massless quark field, h is a heavy-quark field 
defined in HQET, A s is a soft gluon field, and iD£ c = id^ + gA^ c is the covariant derivative 
containing a hard-collinear gluon field. All fields in the above Lagrangian are evaluated at 
point x, except for the soft gluon field in the first term, which is evaluated at x_ = | {n ■ x) n. 
The explicit form of the leading-order Yang-Mills Lagrangian can be found in . 
The terms up to second order in the expansion in %/A are 

/•(i) _ /•(!) , /*(!) , /-CO 

''-SCET — Mj MfM ' 

M2) _ r {2) r {2) r {2) r {2) (A s 
''-SCET — Mj /i % ' I V 

where 



(2) 



2nib 



h(iD s ) 2 h + ^ha, v gG^h 



(5) 



is the next-to- leading term in the expansion of the HQET Lagrangian jTHj. Expressions for 
the remaining Lagrangian corrections have been presented in [T4*] . 



2 



While it is consistent to apply a perturbative expansion at the hard and hard-collinear 
scales, the soft-gluon couplings to hard-collinear fields are non-perturbative and must be 
treated to all orders in the coupling constant. For instance, the hard-collinear quark propaga- 
tor derived from © should be taken to be the propagator in the background of the soft gluon 
field, summing up arbitrarily many insertions of the field A s . The most convenient way of 
achieving this summation is to decouple the leading-order interactions between the soft gluon 
field and hard-collinear fields in (j3J) with the help of a field redefinition, under which |T3*1 ITS] 

£(*) = 5(x_) £<°>(s) , Al(x) = S(x_) A^(x) &{x-) , (6) 

where 

/ V 

(7) 



S(x) = Pexp ^ig J dt n ■ A s (x + tn)^ 



is a soft Wilson line along the n direction. Introducing the new fields into the Lagrangian 
yields 

4 0) = F> | (*» • *>$ + ¥1 ~. — ^ , (8) 

1 V ln - D hc J 

and similarly all interactions between soft and hard-collinear gluon fields are removed from 
the Yang-Mills Lagrangian £yM- The propagator of the new hard-collinear quark field is now 
given by the simple expression 

At(x-y) = (0\T{Z®(x)P\y)}\0) = ±J ^ ^ c -*<.-«) . 0) 

The effect of soft-gluon attachments is taken into account by factors of the Wilson line S in 
the results below. 

We now list the expressions for the subleading corrections to the SCET Lagrangian in 
terms of the redefined fields, using the formalism of gauge-invariant building blocks [T7j . We 
define 

X = W^ 0) , = W^iD^w) , (10) 

where 



W = Pexp ( ig J dtn- A^(x + tn) J 



(11) 



is a hard-collinear Wilson line. These "calligraphic" fields are invariant under both hard- 
collinear and soft gauge transformations. Note that in the light-cone gauge, n ■ A^. = 0, we 
simply have X = ^ and A^ c = gA^c . In terms of these fields, the results compiled in 
take the form 

•(1) _ <r $ ^ 



4 =X£a>" {S^gG, u S) x X. 
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r (2) 



fi ( n ■ x 



n»n v 



+ X^-(if Xhc - 



,1C I 



7 I (S^gG^S)^ + (S^gG fMU S) x 



in - d 2 



in ■ d 



Lhc 



DC . 



4? = (^J^X + h.C. 



(12) 



where iD^ c = i<9 M +./l^ c , and we have dropped the subscript "s" on the soft covariant derivative 
and field strength. The expression for C^j will not be needed for our analysis. The notation 
(. . .) x _ indicates that, in interactions with hard-collinear fields, soft fields are multipole ex- 
panded and live at position x_, whereas hard-collinear fields are always evaluated at position 
x. Because the SCET Lagrangian is not renormalized |14j . the above expressions are valid to 
all orders in perturbation theory. 

Next, we need the expressions for heavy-light current operators in SCET. In general, a 
QCD current qTb matches onto 



q(x) T b(x) 



-tm^v-x 



+ w + jr + w + j, 



r(i) 



(2) 



K 1 ) 



r(2) 



B 



B 



(13) 



where we distinguish between type- A "two-particle" operators and type-5 "three-particle' : 
operators [TH| . The operators arising at tree level are [T4~l H%1 HU] 



J 



(o) 



XT (S*h) 



in ■ d 



r(^) ; 



(2) 



xr 



n ■ x 



{S^.Dh) x + X -^{S^D,D v h) 



2?m 



and 



J 



(i) 



J 



(2) 
B 



+ X-ii, 



in ■ d 



x^4 ±hc —^r(tfh) x -xr^4 ±hc (s^h), 

z in - d ^ mh 



- XT 

-X 



_l_ + J_)„.^ (stA)i 



^ in - d 



2mi 



(14) 



-xi (s^ + x ^ .\ lrU^h) x _ ■ (is) 



in ■ d 



m b 



in 



d 2 



The expressions for the currents beyond tree level are more complicated, primarily because 
several new Dirac structures appear. The relevant formulae are known at leading and 
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next-to-leading order ^Sj i n the power expansion. The corresponding results for the currents 

(2) (2) 

J A and J B have not yet been derived. They would be needed if the analysis in this paper 
should be extended beyond tree level. 

If perturbative corrections at the hard-collinear scale are neglected, the hard-collinear gluon 
fields can be dropped, and the above expressions for the effective Lagrangians and currents 
simplify. In this approximation X — > A^ c — > 0, iD^ c — > id^, and — > 0. While this 
leads to great simplifications in the calculation, we stress that the structures of soft fields that 
arise do not simplify. We find operators containing S^gG^S, S^[D p , gG^S, S^h, S^D^h, and 
S j; D f j i D u h, and the same operators would arise if the calculation was extended beyond the 
tree approximation. The only exception is that we no longer retain Lagrangian corrections 
containing the soft quark field, because — ► in the limit where the hard-collinear gluon 
field is neglected. (The term qSi^±^ ' is forbidden by momentum conservation.) In Section^ 
we analyze the subleading shape functions introduced at 0(a s ) by two insertions of . 

With all the definitions in place, we are now ready to evaluate the current correlator 
including terms of up to second order in y/X, working at lowest order in a s at the hard and 
hard-collinear scales. The leading term is readily found to be 

T<f> = -J^^-^J ^^±2. (te) r,![r,(stfc)^. (ie) 

First-order corrections in vA vanish by rotational invariance in the transverse plane (provided 
we choose the coordinate system such that v j_ = and q± = 0), i.e., = 0. At second 
order in the expansion the correlator receives several contributions, which can be represented 
symbolically as 

jt(2) J (0) ; Jt(l)j(l) ; Jt(0)j(2) j (17) 



Jtd) J(0) 



d 4 zC { ? 



and 



Jt(0) 



Jt(0) j(0) 



d A zCf 



Jt(°) J(°) / d A z 



^-f I ~h 



d w £ 



(i) 



(19) 



where at tree level only the type-A current operators appear. Examples of these time-ordered 
products are depicted in Figure [TJ At 0(a s ), one must include a tree-level contribution of the 
form 

Jt(0) J(0) 



foe® I tfwcw 



(20) 



shown by the last diagram in the figure. Beyond tree level, one would also have to include the 
contributions from the type-I? current operators. 

Because the external 5-meson states in the definition of the hadronic tensor contain only 
soft constituents, and because the hard-collinear fields have been decoupled from the soft 
fields in the leading-order SCET Lagrangian (JHJ), it is possible to contract all hard-collinear 
fields in the time-ordered products. At tree level, all we need is the hard-collinear quark 
propagator Derivatives acting on hard-collinear fields give powers of p in momentum space, 
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hli 

*___>___0 fay » 

i \ A A /fu 

Figure 1: Representative examples of time-ordered products contributing to the 
current correlator in SCET. Double lines show heavy-quark fields, dashed ones 
hard-collinear fields, and wavy lines denote the external currents. Lagrangian inser- 
tions and higher-order effective current operators are exemplified by soft gluons. 





whereas components of x M appearing in the multipole-expanded expressions for the effective 
Lagrangians and currents can be turned into derivatives d/dp acting on the momentum-space 
amplitudes. Each insertion of a SCET Lagrangian correction introduces an integral over soft 
fields located along the n light-cone. Consequently, the time-ordered products in (fTTj) lead to 
expressions involving bi-local operators as in (fTfij). while those in (JTBJ), (fTUj). and (|2U|l also lead 
to tri- and quadri- local operators. At first sight, this would seem to require the introduction 
of complicated sub leading shape functions depending on up to three momentum variables u>i, 
defined in terms of the Fourier transforms of the matrix elements of the non-local operators. 
However, the non-localities can be reduced using partial-fraction identities for the resulting 
hard-collinear quark propagators. At tree level, it suffices to define shape functions of a single 
variable to. 

To see how this works, consider the effect of an insertion of the Lagrangian Cf in G2D. 
Since all hard-collinear fields must be contracted, we may consider without loss of generality 
the expression 

L( yi ,y 2 ) = ijd 4 z At{ Vl ~z)\zl n v {S^gG, v S) z _ A ( (z - y 2 ) . (21) 

The fact that the soft fields live at position z_ implies that the two hard-collinear quark 
propagators carry the same momentum components n-p and p±. In analogy with the definition 
of the hard-collinear calligraphic gluon field, we now introduce the soft field [TJj 

o 

A slx {x) = S\x) {iD^S){x) = - J dtn v (S^gG^S) (x + tn) , (22) 

— oo 

which allows us to write n u (S^ gG^ v S) z _ as a derivative, — n ■ d z A sfl (z_). Integrating by parts 
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in (|2H), we find 



f $ 

L{yi,y 2 ) = i / d zAt(yt - z) ~z£A Bli (z-)n ■ d z A t {z - y 2 ) 

+ ij d 4 z[n- d z Az{ yi - z)\ | zl A slx {zJ) A^z - y 2 ). 
The hard-collinear propagator is a Green's function obeying the differential equation 



n ■ d + 



3_ 

n ■ d 



AAx) 



( J(4)( 



X 



(23) 



(24) 



This allows us to write 

L(vi,y2) = 



y 2 ± ■ A s {y 2 -) - y-LL • A a {yx-) A^{y x - y 2 ) 

-iftfiz As( yi -z)£z ± - A s {zJ) ^L-A^z-y 2 ) 
J 2 n ■ d z 



i I d z 



d 2 

°z± 

n ■ d. 



A € (2/i - z) 



- z± ■ A s (z-) Az(z - y 2 ) 



(25) 



The terms involving transverse derivatives vanish at tree level, because they provide powers 
of transverse momenta of external lines, which are zero (recall that v ± = and q± = 0). This 
leaves the terms shown in the first line, in which the z integral has been eliminated, and in 
which the product of two propagators has been reduced to a single propagator. In momentum 
space, these manipulations correspond to the partial-fraction identity 



n ■ pn ■ p + p 2 ± n ■ p (n ■ p + u) + p 2 j_ n-puo 



n-pn-p + p 2 ± n ■ p (n • p + oS) + p\ 



(26) 

valid for two momenta that differ only in their n ■ p components. 

In the sum of all terms many cancellations and simplifications take place, and we find the 
rather simple result 



Tjp = - I d 4 xe i{q - m » v) - x 



d P Ap-x n 'P 



(27T) 



p 2 + ie 



(27) 



n=l 



where 



Oi(x-) 

o 2 (x_; 

OAx-, 



j d'zT{{hS) Q T^T 3 {tfh) x Ct\z)}, 

(^) r,|r, (sUph) x _ + (hi-ip^r^r, {tfh). 



2mb 
1 

n ■ p 



n 



n-x/2 



(28) 



In deriving these expressions, we have made use of the identity iD^ = S(id^ + A%)S^ with A% 
as defined in (j22j). Note that on the space of forward matrix elements the result for is 
hermitean, because we can integrate by parts and use translational invariance. 

In the last step, we can simplify the operator O2 by noting that the HQET equation of 
motion, iv ■ Dh = 0, along with i/>h = h, implies 

S ] ilph = S l ip ± h + (i-i)in- d(S ] h) . (29) 

It follows that all gauge- covariant derivatives of the heavy-quark fields are perpendicular 
derivatives. This fact restricts the number of subleading shape functions. 



3 Definition of subleading shape functions 

The Dirac structure of the operators O n can be simplified noting that the heavy-quark fields 
h of HQET are two-component spinor fields, so that between h. . .h the Dirac basis collapses 
to a set of four basis matrices (l,cr). In four-component notation, these are the upper-left 
2x2 blocks of (1,775) (in the Dirac representation). Instead of 775, we are free to take the 
matrices 7^75 and ^75. It follows that, between the P v = ~(1 + ^) projectors supplied by the 
heavy-quark fields, any Dirac matrix T can be decomposed as 

r - ~ tr (r p v ) 1 - i tr [r p v a - i) 75 ] i l5 - i tr (r p v 1±lll5 ) 7^75 • (30) 

We denote by 

(B(vU...h\B<v)) (31) 

the forward S-meson matrix element of any HQET operator. Rotational invariance in the 
transverse plane implies that transverse indices can only be contracted using the symmetric 
and anti-symmetric tensors (we set eoi23 = 1) 

_|_ n ^ r) ^ 

flf = fT- — 2 , <% = er#v a np. (32) 
It follows that the only non- vanishing matrix elements are 

((hS) (S^h) x _) = J 
(z / d 4 zT{(hS) (S^h) x Cf\z)}) = — [ d^e-H*^), 

J TTlff J 

■ pp. i> 

((hS) rtYB(&iDlh) x _) = -\\ dwe-Wt( w ), 
dt({hS) Q {SUD^DlS) tn {&h) a _) = 9 -f J due-^u(u;), 



n-x/2 

—i 



o 

n-x/2 



dt({hS)J l5 {SUD^D» ± S) tn {tfh) x _) = J due-^™v(u). (33) 



o 
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If radiative corrections at the hard scale are included, it would be more appropriate to split 
up s(u) = Skm(^) + Cma g s mag (uj), where C mag is the Wilson coefficient of the chromo-magnetic 
operator in the subleading HQET Lagrangian. This ensures that the shape functions remain 
independent of the heavy-quark mass. The definitions of the functions t, u, v are chosen such 
that 



(ft(0)^[0,ar_] (ip±h)(x-)) = J due-^ n - x t(u), 

n-x/2 

-i [ dt(h(0) [0,tn] (iD ± ) 2 (tn) [tn, X -]h{xJ)) = [ dcu e~^ x u(u) , 



—i 



o 

n-x/2 

ifl 



dt (h(0) L [o, in] a^ v gG^(tn) [tn, z_] h(x.)) = \ duo v(u) , (34) 



o 



where [x,y] = S(x)S' f (y) is a product of two infinite-length soft Wilson lines, which on the 
light cone (i.e., for x,y \\ n) collapses to a straight Wilson line of finite length connecting x 
and y. 

We also need a variation of the first matrix element in (J34|) . in which the derivative is 
located at position 0. Using hermitean conjugation, translational invariance, and the reality 
of t(u), which follows from parity and time-reversal invariance of the strong interactions, we 
find that 

((hifi±)(p)1i[0,xJ\h(x-)) = {h(p)t[Q,x-](ip±h)(x-)) , (35) 

implying that all terms containing a single insertion of D± can be related to the function t(u). 
From this relation, it follows that 

due-^ n -*t{u) = ((hS) Q £ [4s±(x-) -4s±(0)] (^h) x _) 

n-x/2 

■A 

dt (h{0) | [0, tn] ^n u gG av {tn) [tn, x_] h(xJ)) , (36) 

which defines the function t(u) in terms of a matrix element of the field-strength tensor. 

It is now straightforward to express the forward matrix element of the current correlator 
Tij in terms of shape functions. The resulting traces of Dirac matrices can be simplified using 
identities for the tensor derived in [20] • Taking the imaginary part, we obtain for the 
hadronic tensor 

Wif = [ duj 5{n ■ p + uj) S(uj) Ti , WW = , (37) 



= I dcu5(n-p + cu) 



uS(u) + t(u) s(u) t(uj) u(uj) v(w) ' 

±2 H J-1+ J- 3 + — J-l — — 1 4 

rrib nib n ■ p n ■ p n ■ p 
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where now p = m^v — q, and 



To 



1 



tr 



1 

-tr 



r, t~ 



3 2 



tr 



r J- r p 
T i lo 7s Ii — ^— 7^75 



T 4 = - tr 



r f ^75 



J 2 
1 + V 



- ?i) 75 



(3? 



It follows from (|37|) that the subleading shape functions s(u;) and u(uS) always come with the 
same trace as the leading shape function S(u). However, u(u) is divided by the kinematic 
variable n ■ p, and so it does not enter in a universal (i.e., process- independent) combination 
with S(u). 

We can now specialize our result to the case of semileptonic decay, for which Ti = 7^(1—75) 
and Tj = 7^(1 — 75). This yields 



dcu S(n ■ p + to) < (n"v v + n v v" - g" v - %^ va<i n a v p ) 



nib ) rrib n ■ p 



2(n'V + n v v") ^ + 2n»n v 



n ■ p 



uj S(cu) 



t(u) t(u) +v(uj) 



n ■ p 



(39) 



Similarly, for the contribution of the dipole operator Q-j 1 to B — > X s 7 decay, the Dirac 
structures are (up to prefactors) = |[7j\$](l — 7s) an d — z[$>7i.](l +7s)j where the 
indices /i, v are contracted with the transverse polarization vector of the photon. In this case 
we obtain 



W v = ^ _ g ^ duo5(n-p + oj) 



1 _ s(u) + s (^) - *M + u(u)-v(u) 
nib J rnb n ■ p 



(40) 



We stress, however, that while at leading power in AQ CD /m fe the dipole operator gives the 
only tree-level contribution to the B — > X s ^ decay rate, this is no longer the case when 
power corrections are included. For instance, interference terms of the dipole operator with 
current-current operators can lead to new subleading shape functions even at lowest order in 
perturbation theory. To derive these structures, it would be necessary to match the entire 
effective weak Hamiltonian for B —>■ X s ^ decay onto SCET operators [6]. This task still has 
to be completed beyond the leading order in A. Contrary to claims in [HI EI], a complete 
description of tree-level subleading shape-function effects in B — > X s ^ decay is therefore still 
lacking. 



4 Moment relations and comparison with the literature 

Moments of the shape functions can be related to forward 5-meson matrix elements of lo- 
cal HQET operators 0. In particular, setting x = in the defining relations (|3*3*|) yields 
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expressions for the normalization integrals of the shape functions. They are 



dui S(uS) 



dui {s(ui),t(ui),u(ui),v(ui)} = . 



(41) 



The vanishing of the norm of all subleading shape functions is a consequence of Luke's theorem 
|2~T| . and it ensures that there are no first-order AQ CD /m b corrections to total inclusive decay 
rates. For the functions t, u, v this is an obvious consequence of the fact that the integration 
domain in ()36|) and (}3"4*j) shrinks to zero in the limit x — > 0. The interpretation of ()41|) is that 
subleading shape functions lead to local distortions of inclusive spectra, which cancel out when 
the spectra are integrated over a sufficiently large region in phase space. The first moments 
characterize the strength of the distortions, while higher moments determine their shape. 

Taking a derivative in ■ d x in the definitions (|3*3|) brings down a factor of ui under the 
integrals on the right-hand side. Setting then x — > yields a set of relations for the first 
moments of the shape functions. The resulting matrix elements can be evaluated by means of 
the relations 



and |231 



{hY aP iD a iD p h) = itr 



(i / <PzT{(hiD*h)(0)C%>(z)}) = -^<>C(0)> = -v 



(9 



Ai 



)2L + ia ^h 

I Q O 



1 + 



(42) 



(2), 



M Ai + 3C mag A2 



2mi 



(43) 



where Ai and A2 are the familiar HQET parameters arising in the parameterization of second- 
order power corrections to inclusive decay spectra, and C mag = 1 at tree level. We obtain 



(kou s{u) = — 



A x + 3A S 



2Ai 



dujujt{uj) = —X 2 



dujuv(uj) = A 



2 • 



(44) 



We also recall that the first two moments of the leading shape function are J duiuj S(uj) = 
and / duuo 2 S{uj) = -Ai/3 0. 

It is appropriate at this point to relate our subleading shape functions to those defined 
by Bauer et al. jHl CHj. Using the notations of their second paper, we find S(u) = /(— uj) at 
leading order, and 



s(u) 



-uj 



t(u) 



-/ill 



-UJ) 



U(UJ) 



-Go 



-UJ) 



VlUJ 



-Ho(-Lo) 



(45) 



for the subleading functions. Their first paper uses a dimensionless momentum variable, so 
the relations are S(uj) = f(—uj/mb)/mb etc. (In comparing, one must take note of the fact 
that we use a different sign convention for the Levi-Civita tensor. Also, we believe there are 
typographic errors in the definitions of the functions g 2 and h 2 in [S].) 
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Comparing our result (J37|) for the hadronic tensor with the corresponding expressions ob- 
tained in [SJ [""0] we find some sources of disagreement. As mentioned earlier, the function 
u(lu) is divided by n ■ p, whereas Bauer et al. state that it enters in the universal combination 
S(u) + u(u)/mb with the leading shape function. Also, it appears to us that the term pro- 
portional to the trace T 3 , which adds a contribution to the last structure function in ("""*"]) . is 
missed by these authors. 

5 Contributions from four-quark operators 

The last diagram in Figure ["J shows a tree-level contribution to the hadronic tensor involving 
two insertions of the subleading SCET Lagrangian . The exchange of a hard-collinear 
gluon implies that this graph is of order g 2 , and so it vanishes in the limit where a s is set 
to zero. The fact that the suppression factor is 7ra s instead of ot s /ir reflects the phase-space 
enhancement of four-quark tree- level graphs compared with loop diagrams [21]. One might 
therefore expect that the four-quark contribution is numerically as important as the other 
tree-level subleading shape-function contributions. 

Applying the partial-fraction identity (|26|) twice, we find that the resulting contribution 
to the correlator T-^ can be written as in ("2T]) . adding a fifth operator to the sum. It reads 

n-x/2 n-x/2 

°^ x -) = iri I dh I dt 2 {hs)j t i^t a (si q ) ti jqs) t2n y ± ir j t a (tfh) x , (46) 

n ■ p j j 

ti 

where t a are the generators of color SU(iV c ). Note that the field insertions are ordered according 
to "light-cone time" x_, just as they appear in the Feynman diagram in Figure ["3 This is 
a general result. Because the minus components n ■ ph c of hard-collinear momenta are large, 
of order m^, hard-collinear fields always propagate forward in light-cone time. Turning, for 
instance, a forward-moving hard-collinear quark into a backward-moving hard-collinear anti- 
quark would require a hard quantum fluctuation, which is already integrated out in SCET. 
As a result, Feynman amplitudes in SCET are ordered with respect to light-cone time, and 
that ordering is preserved in the matching onto HQET. This discussion explains why all our 
operators have the property that the coordinates z_ of soft fields range from to x_ in an 
ordered fashion. The results can therefore always be expressed in terms of bi-local operators 
depending only on and x_. However, at present we cannot exclude the possibility of non- 
trivial weight functions under these integrals, which could arise at higher orders in perturbation 
theory. If present, they may require a generalization of our definitions of subleading shape 
functions. 

Returning to the case of the four-quark operator in (J46|) . we note that its contribution 
vanishes in the vacuum-insertion approximation due to the color-octet structure of the heavy- 
light quark bilinears. While this approximation is admittedly naive, phenomenological ev- 
idence based on studies of 5-meson lifetimes |24[ 12*5] and theoretical arguments based on 
lattice calculations "2UJ [2Z] and QCD sum rules [2H] support the notion that matrix elements 
which vanish in the vacuum- insert ion approximation are numerically suppressed, typically by 
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an order of magnitude. It is thus very unlikely that the four-quark operator in (J46|) could give 
a larger contribution than loop-suppressed 0(a s ) corrections, which we have neglected. Note 
also that the contribution of the corresponding local operator to total inclusive rates, which is 
what remains when the corresponding shape functions are integrated over a sufficiently large 
domain, is bound to be tiny. The effect is of order (Aqcb / m>b) 3 , and it is a s -suppressed with 
respect to a four-quark contribution discussed by Voloshin [29 3 , whose effect on the total decay 
rate is believed to be at most 3%. 

Nevertheless, it is interesting to study the structure of the operator 5 (x_) in more detail 
and define corresponding subleading shape functions. The decomposition (|3T)j) implies that 
the Dirac structure can be rearranged in the form (omitting factors of S, and color indices 
for simplicity) 



\ 



ii5 



\ 



(1/1 - ^75 
V 7^75 / 



rv 



lp H 1 



(47) 



where the notation implies a sum over the three rows in the equation. In the next step, we 
use that between ifi . . .ifi any Dirac matrix can be decomposed as 



r - \ tr(r tr(r i l5 ) \ 75 - \ tr(r i^) 1 7^ 



(48) 



Only the first two terms survive after contraction of the index p in (|47|) . At this stage we are 
left with four-quark operators with Dirac structures 



hhq^(^ 5 )q, hifi^ h qiji (75) q , h 7^75 h q $ (75) q ■ 



(49) 



where (75) means 1 or 75. Lorentz and parity invariance imply that only the two Lorentz-scalar 
structures hhqifiq and h ^75 h q ^75 q can have non-zero forward matrix elements between B- 
meson states. Putting back color factors and soft Wilson lines, we define the corresponding 
subleading shape functions 



2(-i) 



n-x/2 n-x/2 

d tl J ^([(^) tJJta(^) :c _] z [(g5) t2n ]^[(^) ti J fc ) 

fi 



n-x/2 n-x/2 

2H) 2 / dt x J ^ 2 ([(^) t a ]^7 5 [t a (^) x JJ(g5) t2 J^7 5 [(^) ti J^ 
ti 

dwe-i^ x f v (u), 



(50) 
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where k, I are color indices. Using the same notation as in (|37|) . the corresponding contributions 
to the hadronic tensor are given by 



(2) 



-na. 



■lq 



duj 8{n ■ p + uj) 



1 1 + — i4 



n • p 



(51) 



with the same traces as defined in (J38J). As far as the spin structure is concerned, the four- 
quark contributions can thus be absorbed into a redefinition of the subleading shape functions 
u and v, namely 



u(u>) = u(uj) - ira s f u (u) , v{uj) = v{uj) + 7ca s f v (uj) 



(52) 



Note that the definitions (JH0|) imply that the normalization integrals as well as the first mo- 
ments of the functions f u (u) and f v (u) vanish (because the integration domain is of second 
order in n ■ x), and therefore the new functions u and v have the same first moments as the 
original ones, see (j44j). 

While the shape functions S, s, t, u, v are expected to be identical for charged and neutral 
B mesons up to tiny isospin-breaking corrections, this is no longer the case for the four-quark 
shape functions f u and f v . The values of these functions will depend crucially on whether the 
light-quark flavor q in the four-quark operators matches that of the 5-meson spectator quark. 
In the semileptonic decay B — ► X u l~u, the difference between the subleading shape functions 
f u and /„ for B~ and B° mesons is likely to be one of the dominant sources of isospin-breaking 
effects on the decay distributions. 



6 Applications 

We are now ready to study the phenomenological implications of our results. We absorb 
the contributions from four-quark operators into the functions u(u) and v(u) defined in (JH2J). 
The moment constraints on the shape functions can be summarized by their expansions in 
distributions [1], which read 

S(cj) = 6(u>) - ^ S"(u) + ..., s(u) = Al+ 3A2 6'(u) + ..., (53) 
o 2 

2A 

t(u) = X 2 5'(co) + ... , u(uj) = — -Ls'( u ) + ..., v(u>) = -X 2 S'(co) + ... . 

These expressions allow us to test our results against existing predictions for inclusive spectra 
obtained using a conventional heavy-quark expansion. 

The analytic properties of the shape functions are such that they have support for — oo < 
uj < Aqo, where A ra is the asymptotic value of the mass difference {m B — m b ) mb ^ 00 in the 
heavy-quark limit. This parameter differs from the physical value of A by power-suppressed 
terms, 

— — Al ~\~ 3A2 fr a\ 

A = m B -m b = A 00 h . . . . (54) 

Zrrib 

We would like the support in u to extend over the physical interval — 00 < uj < A, since this 
will ensure that the kinematic boundaries for decay distributions take their physical values 
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set by the true 5-meson mass. This can be achieved by shifting the arguments of all shape 
functions by a small amount Au = ~(Ai + 3\ 2 )/m b . For the subleading shape functions this 
changes nothing to the order we are working, since t{u + Au) = t{u) + . . . etc., where the 
dots represent terms of higher order in l/m b . For the leading-order shape function, however, 
this shift produces a new l/m b correction: S(u) = S(u + Au) — Auj S'(u + Au) + . . . , where 
the prime denotes a derivative with respect to the argument. Using the fact that S(u) and 
s(u)/m b always appear together, we can absorb the extra term into a redefinition of the 
subleading-shape function s, defining a new function 



s (u)=s(u)--(X 1 + 3X 2 )S'(u) 



(55) 



From (gH) , 
definitions, 



it follows that the first moment of the function sq vanishes. In terms of these 



S(u + Acu) = S{u + Au) + 



sq(uj + Au) 



s(u 



S(u) + + 



(56) 



m b m b 

In other words, our expressions for the hadronic tensors in and remain valid when 
in all shape functions the argument is shifted from u to u + Au, except that the subleading 
shape function s must be replaced with the redefined function sq, whose norm and first moment 
vanish. Once this is done, the integrals over u extend from — oo up to the physical value of A. 

It is now convenient to introduce a new variable P + = n • p + A = mg —n - q = Ex — \Px \ , 
which is the plus component of the total momentum of the final-state hadronic jet, and to 
express the shape functions as functions of u = A— u 0. They have support in this variable for 
< u < oo. The 5-functions in the tree-level expressions ([39)1 and (|40|) for the hadronic tensors 
set u = P + . We denote functions of u by a hat, e.g. S(u) = S(A — u + Au) = ^(Aoo — u), 
t(u) = t(Aoo — u), and similarly for the other functions. In the equations below, A always 
refers to the physical parameter defined with the true 5-meson mass. 

We start by presenting the triple differential rate for the inclusive semileptonic decay 
B — > X u l~v in terms of the kinematic variables 

2Ei 



P + , y = ^ = l-^, x = l-^i. (57) 
m,b m b m b 

Using the expression for the hadronic tensor in (|39~j) combined with the rate equations derived 
30 , we obtain 



n ■ q 



in 



1 



ci 3 r 



r dx dy dP+ 



12(3/ -x) {(1-y + x) 



1 + 



2(A-P + 



5(P + 



i(p+) , u(p+) - «(p 4 



m b 



ym b 



+ 1 



2x\ A - P. 



m b 



S{P + ) 



2x{\- 2y) t{P + ) 2x v{P. 



+ 



m b 




(58) 



This expression is exact at tree level and to order Aqcj)/m b in the heavy-quark expansion. 
The phase space for these variables is such that 



< P + < m B - 2Ei = m b x + A . 



m b 



< x < y < 1 



(59) 
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and the hadronic tensor has support only for y > 0. Above we assume that x and y are of 
0(1). The collinear expansion breaks down in the region of small y. A restriction to small 
x is allowed and corresponds to the so-called "endpoint region" of the charged-lepton energy- 
spectrum. In this case, however, it is necessary for consistency to expand the expression above 
in powers of x = O '(Aqcd / 'mb) ■ 

As a first application, we integrate ()58|) over the lepton energy, finding 



1 d 2 T 



T dy dP + 



2y 2 {3-2y) 



+ 22/(6 - 5y) 



1 + 



2(A-P + 

m b 



S(P+) + 



i(P+) , u{P. 



+ 



ym b 



^±S(P + )+2y(l-2y) 
m b 



2t(P + )-v(P + 
m b 



(60) 



The variables y and P + allow us to reconstruct any kinematic quantity characterizing the 
final-state hadronic system, such as the total hadronic energy, 2E X = P+ + ym b + A, or the 
invariant hadronic mass, M\ = P + (ym b + A). Of particular interest for a measurement of 
\V ub \ are the spectra in the variables P + and Mx- For the distribution in P + , we obtain 



1 dT 



r dp. 



14(A - P + 
3m;, 



S(P+) + 



t(P + )+5u(P + )+v(P+ 
3m b 



(61) 



Next, denoting sh = M x and A s = sh/ttib, we find for the hadronic invariant mass spectrum 



1 dT 



T ds 



H 



1 

m B 



dP+ 



/(P + ,A S /P + 



(62) 



where 



f(P + ,r) = 2r 2 (3-2r) 



1 + 



- 2r 



2(3A - P + 
m b 

A 



S(P+ 



t(P + . 
m b 



+ 



u(P+ 



(6 - 5r) — + 2r 

m b m b 



S{P+) + 2r(l - 2r) 



rm b 

2i(p + ) - 



0(P 4 



m b 



(63) 



A comment is in order concerning the fact that the integration over P + in ()62)1 is extended 
to infinity, while kinematically P + < The point is that near its kinematic limit the 

momentum P + ~ a/^b^V^cd is no longer of order Aqcd ; and hence the collinear expansion 
breaks down. However, the contribution to the rate resulting from the vicinity of this region 
is suppressed by more than a single power of A QCD /m b jl] . It is therefore required to set the 
upper limit to infinity so as to avoid spurious higher-order power corrections. Once this is 
done, it is theoretically consistent to treat P + as a quantity of order Aq CD . 

As a final application, we integrate (JSHJ) over the hadronic variables y and P + to obtain 
the charged-lepton energy distribution 



1 dT 
f~dE, 



du 



m B 



1 + 



6 (A 



m b 



i(u>) + 3m(o>) — 3v(u) 



m b 



(64) 
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For purposes of comparison with previous authors, we also give the result for the dipole- 
operator contribution to the B — > X s ^ photon spectrum, stressing once again that this does 
not provide a complete description of all tree-level subleading shape-function effects. Taking 
into account that the photon spectrum is proportional to the expression in (j4*Uj) times a factor 



E^, we find 



1 dr 



2T dE^ 



1 + 



2(A-A 

m b 



-t(P + )+u(P + )-v(P+ 
m b 



+ ... , 



(65) 



where P+ 



2E~ in this case. 



m B ^ 7 

Despite the discrepancies mentioned earlier, we agree with [H] on the result for the B — > X s 7 
photon spectrum (having, however, added the new contribution from four-quark operators), 
and we also find complete agreement with a recent calculation of the hadronic invariant mass 
spectrum for semileptonic B — > X u l~v decay by Luke et al. However, we do not confirm 
the result for the charged-lepton energy spectrum obtained in JU]) when converted to our 
notation. These authors have (apart from an overall missing factor of 2) a coefficient of 2 
instead of our 6 in front of the term proportional to (A — Co), and a coefficient of 1 instead of 
our 3 in front of u(Cj). Our results can be tested against established formulae by using the 
moment expansions in We find that (JBTj). (JEI|, and are consistent with expressions 

for the corresponding spectra derived in [7j, j2H|, and jT], respectively. Amusingly, the result 
for the lepton spectrum found in |IUj is also consistent with the moment expansion. The 
reason is that this spectrum differs from our result by a term involving the combination 
4(a;/m&) 5" (a;) + 2u(uj), in which the 5' (ou) terms from (Jo"H|) cancel each other. 

Of particular phenomenological importance are various fractions of all B — > X u V~v decays 
that pass certain experimental cuts, which are chosen so as to eliminate (or reduce) the 
background from semileptonic decays with charm hadrons in the final state. The most common 
cuts are E t > E with A E = m B — 2E < m 2 D /mB ("lepton endpoint"), sh < s with A M = 
s o/ m B < ^d/itt^b ("low hadronic mass"), and P + < Ap with Ap < m 2 D /mB ("low hadronic 
plus momentum"). We define corresponding fractions Fg^A^), F M (A M ), and F P (A P ). While 
the result for Fp(Ap) is simply an integral of the spectrum in (JBTj) over the range < P + < Ap, 
the expressions for the other two event fractions are more complicated, and we quote them 
here for completeness. We obtain 



Fe{Ai 



du 



2{A 1 



nip — uj 



6(A -Cj)\ A . _ i(u) + 3u(tu) - 3v(u) 
b{u) H 



m b 



(66) 



and 



Em (Am) = J dP+ 
o 

oo 



I + 14(A-P + ) ^ g + t(P + ) + 5u(P+) + v(P+ 



3m b 



+ J dP + g(P + ,A M /P + ), 



3m b 



(67) 
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where A M = s /m B , and 



g(P+,r) =r 3 (2-r) 



2(3A-P+)\ - i(R 



1 + ^ S(P + ) + 



+ > 



m b l m b 



2r< 

T 



9 - 5r — + 2r — 



+ ^ (9 _ 4r )^) + r! ( 3-4r) 2 ^-^. (68) 
3 m& 3 mi 

In order to illustrate the numerical impact of subleading shape functions, we focus on the 
single differential spectra in (|61|h ([62)1 . and in semileptonic I? decay. For the leading 
shape function we use the parametrization 

^ode 1 (^) = T |y^e-^/ A (69) 

with A = 0.63 GeV and b = 2.93, corresponding to the default choice in We employ 
two models for the subleading shape functions based on the moment expansion in (|53)L one 
where S'(u)) is replaced by a derivative — S' model (u) of the leading shape function, and one 
where it is replaced by [3(A — u))/Ai] S'modei(^)- Both models satisfy the moment constraints 
identically, while giving rather different shapes for the subleading structure functions. We take 
A2 = 0.12 GeV 2 , m b = 4.65 GeV, and determine Ai = —0.41 GeV 2 from the second moment of 
the model function (JB9*|). 

Results for the three spectra and the corresponding integrated event fractions are presented 
in Figure 121 In each plot, the dashed-dotted line shows the leading term in the heavy-quark 
expansion, the dashed gray curve includes the effect of the "kinematic" power corrections pro- 
portional to S(u), which are predicted model independently, and the two solid lines include 
the model predictions for the subleading shape-function effects (the lighter of the two lines 
corresponds to the first model). While the effects are clearly visible on the spectra themselves, 
the integrated rate fractions with a cut on P + or sh receive only moderate power correc- 
tions, which alter the leading-power results by typically 10% or less in the region near the 
charm threshold (located at m 2 D /mB ~ 0.66 GeV). Note that the kinematic corrections alone 
overestimate the effect, especially for the case of the P + spectrum. 

The situation is rather different for the charged-lepton energy spectrum, where power 
corrections can have a drastic effect and tend to reduce the event fraction by a large amount. 
This effect can be traced to the particular weight function under the integral in (J6I3J) [TT] . 
Even for large values of 00, the effect of subleading shape functions does not diminish, because 
the linear term in the weight function picks out the first moments of the subleading shape 
functions. As a result, there is a significant theoretical uncertainty in the prediction for the 
event fraction Fe{Ae) near the charm threshold, which imposes limitations on the precision 
with which \V ub \ can be determined from the charged-lepton endpoint region. 

No such large effects are seen for the event fractions with a cut on P + or hadronic invariant 
mass. With Ap t M ~ m 2 D /mB, about 65-85% of all B — > X u l~v events are retained with these 
cuts, and the corresponding fractions are enhanced further when radiative corrections are 
included Once the efficiency is so high, possible uncertainties due to four-quark operator 
contributions ("weak annihilation") become negligible. Our analysis supports the study of 
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Figure 2: Predictions for the P + , sh, and E\ spectra (upper plots) and event frac- 
tions (lower plots). The dash-dotted curves refer to the heavy-quark limit, while the 
solid lines include subleading shape-function effects. The gray dashed curves include 
only those power corrections proportional to the leading-order shape function. 

theoretical uncertainties presented in [7] , where the uncertainty due to power corrections in the 
value of Fp(Ap) near the charm threshold was estimated as 10%. It provides further credibility 
to the strategy of obtaining highly efficient, precise measurements of \V u b\ using experimental 
cuts on P + or sh- We recall at this point that the P + spectrum offers the additional advantage 
that it can be directly related to the B — > X s ^ photon spectrum, cf. (JoTj) and whereas a 
prediction of the hadronic invariant mass spectrum requires knowledge of the shape function 
over a range not accessible in radiative B decay. Note also that the hadronic invariant mass 
distribution receives fractional power corrections of order (Aqcd/^6) 7 with 7 « 1.35 from a 
region of phase space where the collinear expansion breaks down 0]. While these effects are 
formally beyond the accuracy of the present work, they introduce an additional theoretical 
uncertainty that is difficult to quantify (unless they can be shown to cancel). 

The smallness of subleading shape-function effects for some observables can be understood 
based on the moment expansions in (|53j) . For instance, in the P + spectrum (|61|) the first 
moments of the functions t and v cancel in the combination t + v, and the first moment of 
|u cancels against that of y (A — P+) S, leaving a small residual term — g XiS'(P + — A). 
In weighted integrals over the P + spectrum, this can give rise to small second-order power 
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corrections. Higher moments may spoil these cancellations, but their effects are only of order 
(Aqcd / 'mj,) 3 or higher when the spectrum is integrated over a sufficiently large domain. Similar 
cancellations occur for the hadronic invariant mass spectrum (|62j) . and for the photon spectrum 
(|65|h for which the first moments cancel entirely. In all three cases, we have checked that 
taking a mix of the two models for the subleading shape functions, in which different higher 
moments prevent perfect cancellations, does not change the results by a significant amount. 
The case of the charged- lepton spectrum (J64j) is again different. The first moment of 3m cancels 
against that of 6(A — Cj) S. However, the sum t — 3v has a rather large first moment given 
by — 4\ 2 5'(u — A), which in conjunction with the linear weight factor in (J55j) gives rise to a 
significant power correction to the event fraction F E EH E] . 

7 Conclusions 

With the development of soft-collinear effective theory (SCET), theoretical predictions for 
inclusive S-meson decay distributions have recently received renewed attention. While these 
processes had been studied extensively during the past decade, the ever increasing precision 
of the measurements at the B factories requires a new level of accuracy. This is an area where 
SCET is not only a useful conceptual tool, but where it provides concrete means of pushing 
the limits of theoretical calculations. 

Earlier this year, two groups have performed systematic analyses of short-distance correc- 
tions in the semileptonic decay B — > X u l~ v jUIH] and radiative decay B — > X s ^ [6], including 
a complete scale separation and resummation of Sudakov logarithms. As a result, the behavior 
of the leading shape function under renormalization is now well understood. This leaves power 
corrections to the heavy-quark limit as the principal source of theoretical uncertainties. 

In the present work, we have used the formalism of SCET to perform a systematic study 
of such power-suppressed effects. At tree level, the results can be expressed in terms of a set 
of subleading shape functions defined via the Fourier transforms of forward matrix elements 
of bi-local light-cone operators in heavy-quark effective theory. We have identified a new 
contribution arising from four-quark operators, which was not considered previously. We 
have shown that, when shape functions appearing in process- independent combinations are 
combined into single functions, then a total of three subleading shape functions are required 
to describe arbitrary current-induced decay distributions of B mesons into light final-state 
particles. While subleading shape-function effects had been studied in the past, our results do 
not agree with those found in the original papers [HI HU] • 

In the last part of this work, we have presented analytical expressions for a variety of 
distributions in B — > X u l~v decay, which can be used directly for the analysis of experimental 
data. We have also given a formula for the triple differential rate, which allows for arbitrary 
cuts on kinematic variables. While this concludes the problem of tree-level power corrections in 
semileptonic decay, we have stressed that no complete (tree-level) analysis of power-suppressed 
corrections to the B — > decay exists to date. The formalism developed in this work can, 
however, readily be extended to this case. 
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